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Abstract—The finite element method, using flat shell elements, is applied for determining the stress levels induced 
in a typical steam turbine casing due to its self weight. Computer programs have been developed for the analysis 
and frontal solution technique with diagonal pivoting has been used for the solution of the finite element equations. 
Various loading histories corresponding to different orientations of the casing during mechanical handling were 
considered in the analysis. The results of the analysis were compared with limited amount of experimental data and 


were found to be in good agreement. 


INTRODUCTION 


Generally, a medium sized thermal power plant, has a high 
pressure turbine, an intermediate pressure turbine and a 
low pressure turbine. In each turbine, the steam expands 
and the thermal energy is converted to mechanical 
energy for driving the generators to produce electricity. 
The flow diagram of a typical plant is shown in Fig. 1. 
Each turbine has a rotor wherein the deflector nozzles, 
diaphragms, bearings, etc. are housed. In case of the low 
pressure turbine, the casing is a very large structure, 
since it has to support heavy rotor assemblies and to 
handle voluminous amount of steam. A typical LP casing 
for a medium range plant weighs around 30 tons. The 
bottom half of the casing is shown in Fig. 2 and weighs 
around 20 tons. It consists of stiffeners, side plates, 
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conical shell to house the rotor and lifting lugs. The 
various parts of the casing are fabricated separately and 
are welded together in a particular sequence. This 
assembled casing is transported from the final assembly 
shop floor to the site. During lifting of this structure, if 
the proper. instructions are not adhered to, such heavy 
structures are prone to permanent deformations due to 
high stress caused by its own weight if the stiffeners are 
inadequate. The permanent deformation would cause 
difficulty in alignment during erection. 

The purpose of this work is to validate the design of 
one such casing due to the stresses induced during 
handling and to verify the positioning of the stiffeners. 
Since the structure is complex in shape with stiffeners 
connecting the side plates, it is difficult to analyze by 
conventional method. The Finite Element Method (FEM) 
has been used for the analysis. The flat shell elements for 
the side plates and conical portion of the casing and 
space frame elements for stiffeners are used in finite 


Fig. 1. Flow diagram of a thermal plant. 
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Fig. 2. LP cylinder casing. 


element modeling. Three types of loading during lifting 
are considered for the analysis. The FEM analysis is 
checked with the available experimental results. 


FINITE ELEMENT STRESS ANALYSIS 


The LP cylinder casing consists of conical shell and 
flat plates and bars. Hence, it is considered as a shell 
structure and is idealized by an assembly of thin flat 
triangular elements which are assumed to undergo both 
membrane and bending deformations. The stiffeners are 
considered as space frame elements with six degrees of 
freedom at each node. 


Element characteristics 

The triangular flat element, subject to inplane and 
bending forces are shown in Fig. 3. It is assumed that the 
inplane forces do not affect the bending deformation and 
vice versa[1]. Hence, the inplane stiffness matrix [K],° 
and bending stiffness matrix [K],° are assumed to be 
uncoupled. In deriving the element stiffness matrices, 
linear variation is assumed for the inplane displacements 
(u,v) and cubic variation is assumed for the normal 
displacement w. At each node, the displacements u, v, w 
and rotations 6,, 6, and 6, are taken as unknowns. 
Area coordinates are used defining the shape function 
[N] as shown in Fig. 4. The elemental governing stiffness 
equations are obtained by adding the individual stiffness 
matrices due to inplane and bending forces. 

Elemental stiffness matrixes [K],° and [K]}.°. The in- 
plane stiffness matrix [K],° is 6X6 order for the dis- 
placements u, v expressed in terms of the nodal dis- 


placements as 
u = Nii + Niu + Ninn 


(1) 
v = Niv, + Nj t+ Nmtm 

in which 

_ atbix+cy 

NU 9A 

Gi = Xj¥m — XmYj 

bi = ¥i~ Ym 

Ci = Xm — Xj 
and 


A = area of the triangle. 


The stiffness matrix [K],° is computed from ([1], p. 101) 


[K],°= J i [B]"[D],[B] dx dyt 
=(B]"(D],[B]At (2) 
wherein (B] is the strain displacement matrix and [D] is 


the elasticity matrix connecting the stresses {0} to the 
strain {€}, i.e. 


Ox lv 0 x 

Oy S= el v 1 0 €y (3) 
1 

Txy 0 0 oan Yxy 


{a} = [D]p{e} (3a) 
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Fig. 3. Area coordinate system. 


and The bending stiffness matrix (K],° is 9x9 size and is 
obtained from displacement function. 
dul ax 4 ; 
dv ay ‘ w= oaiL1+ a2L2+ a3L3+ aL Lot 3LiL2L3) 
{e}= du, av [> [BH5}, (3b) 2 1 2 1 
a ar + as(Li'L3+3Lh.L2Ls) + ae(L2L3 + 3LiL2L3) 
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+ ao(L3°L2 + Li LoL). (5) 
with 4 
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terms of the nodal displacements are expressed as 


{e}=(QHa} 
= (QUCT {5}.° = [BS}.. (7) 


The generalized “‘stresses” {o} =[Mx, My, Mxy]" and 
strains are related by the elasticity matrix [D], as 


{o} = [D]e fe}. (8) 


The elemental stiffness matrix K,° is expressed as ((1], p. 
236) 


Kp =tcr" (ff te" wwiaieray)icr'. ©) 


Numerical integration scheme suggested by Cowper [2] 
has been used to evaluate the matrix [K].°. 

Assembly procedure. The elemental stiffness charac- 
teristics are evaluated for each element in its element 
coordinate system. 

The element stiffness equations are of the form 


{f}° = [K]*{a} (10) 
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with element stiffness matrix 


Kr. }000:0 
2x2 100 0; 0 
00: Ky $0 
= ie 10 
[K]-s 00: 3x3 {0 (10a) 
wae e---------------! ann 
00:000:0 


The elemental stiffness matrix can be transformed to 
the global system by usual transformation 


([K]& =(T)"[K]*(T] 


wherein [T] is the transformation matrix. 

In the global system, an additional degree of freedom 
will be present because of the transformation. Com- 
patibility is achieved by inserting null row and column in 
the stiffness matrix for the coefficients corresponding to 
the unknown values as in eqn (10a). 

If the elements are coplanar, as is the case in this 
analysis, the final solution will lead to difficulty because 
of the singularity of the matrix. To circumvent this 
difficulty, a fictitious set of rotation stiffness coefficients 
are added in all the elements as suggested by 
Zienkiewicz[1]. For the triangular element, the fictitious 
stiffness coefficients are defined as 


1-05 -0.57)/ 6. 
= aEth - 05 1 -05 He, 
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Fig. 5. Finite Element Discretization of LP casing. 
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Fig. 6. (a) FEM Discretization of central support plate with conical shell. (6) FEM discretization of side plate I. (c) FEM discretization of 
side plate II. (d) FEM discretization of front plate. (e) FEM discretization of rear plate. 


wherein @ is a coefficient and is taken as 0.03 in the 
analysis. The stiffness matrix of the stiffener element is 
given in Appendix 1. 


Discretization and solution procedure 

Even though the structure is geometrically symmetric 
about the axis of the turbine, the loadings due to various 
parts of the assembly, i.e. the weights of the parts which 
are acting on the conical portion of the casing, are not 
exactly symmetric. Apart from this small unsymmetric 
loading, during lifting of the structure, a little tilt would 
alter the symmetricity. Hence, the complete structure is 


considered for the analysis. The structural discretization 
is shown in Fig. 5. The discretization resulted in 438 
elements with 1428 degrees of freedom. The data pre- 
pared for nodal coordinates, and element connectivity 
are checked by computer plot. Each individual com- 
ponents such as conical shell and side plates are plotted 
and they have been checked to avoid any error in data 
(Fig. 6). 

The structure considered in this analysis is a closed 
one and use of any band solver procedure will be un 
economical with respect to computer memory and the 
solution of assembly equations. In order to efficiently 
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Fig. 7. Loading conditions. 


handle the solution of this problem, the frontal solution 
technique[3] is followed. The frontal solution is ad- 
vantageous since the solution procedure proceeds based 
on elementwise elimination. Hence, this is highly suited 
for the analysis of closed structures and the solution can 
be obtained with minimal active core memory. The solu- 
tion procedure presented in [3] has been modified to 
incorporate the diagonal pivoting during the elimination 
process. Incorporation of diagonal pivoting permits use 
of less core and minimizes the input/output operations. 
Also, after assembly of each element, the nodes which 
are not appearing further in any other element are eli- 
minated resulting in more space for assembly of the 
other elements in the solution procedure. 

The computer programs developed for this analysis 
consist of a main program and subroutines for (i) evalu- 
ating the element inplane and bending stiffnesses, (ii) 
frontal solver and (iii) stress computation. The frontal 
solution routine assembles the element equations and 
solves for the nodal displacements. The element stresses 
are computed using the displacements. The stresses in 


the shell elements due to bending are calculated at the 
integration points in both the surfaces. The principal 
stresses are evaluated using the inplane and bending 
stresses. 

The loadings considered for the analysis are (i) the 
perfect lift, (ii) the casing is tilted by 10° about the 
turbine axis (6 = 10°) and (iii) the casing is tilted by 10° 
about the transverse axis (¢ = 10°) (Fig. 7). 

For each loading, the computer time taken for the 
analysis is 56min. of I/O time with 19.5 min. of CPU 
time in a IBM 370/155 systems with 512K bytes 
memory. 

The results obtained for various orientation of the 
casing during lifting are shown in Table 1. 


EXPERIMENTAL ANALYSIS 


An experiment was carried out in the site on the same 
casing by strain gauging. The strain gauges were put at 
selected points as shown in the Fig. 8. The stresses were 
computed using the strains induced by lifting the actual 
model. Only selected points were considered because of 
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Fig. 8. Position of Strain gauges in LP cylinder. 


the restriction due to the capacity of the instrument. The 
stresses computed were given in Table 1. 


Table 1. Comparison of stresses (Kgf/cm) 


6=o=0 6=10 g=10° 
Position FEM Expl FEM FEM 
R2 6.8 6.17 21.09 19.60 
R3 10.7 10.10 23.00 21.00 
R, 2.6 1.98 441 5.10 
Rs 3.1 3.57 2.41 2.80 
Ry 11.2 10.80 26.20 25.1 
Rs 12.6 15.10 29.60 39.9 


RESULTS AND DISCUSSION 


For normal lifting condition for which the experimen- 
tal results are available, the stresses are evaluated using 
FEM modeling and a comparison has been shown in 
Table 1. The variation of the FEM results with the 
experimental values are 4-30% in some locations. This is 
due to the lumping of the weights of the various parts 
which are coming to the conical portion of the casing to 
house the rotor. However, the results are in good 
agreement with the experimental values which provide 
confidence in using the same modeling for different types 
of loading. For the given stiffener arrangements, the 
stresses induced during normal lifting of the structure are 
within the allowable stress values. 

In actual practice, a small amount of tilting may occur 
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about the central line of the casing during mechanical 
handling. This tilting is simulated by simply changing the 
nodal coordinates. The tilting about turbine and trans- 
verse axes have been considered. It is observed that the 
stresses have reached double the value compared to the 
normal lifting conditions for a 10° tilt (Table 1). The 
stresses for this amount of tilt are within admissible limits 
which validates the stiffener arrangements provided. 
However, if the stiffeners are inadequate, the stresses 
will be very high which will result in permanent defor- 
mation at certain locations of the casing leading to 
difficulty in alignment during erection. 


CONCLUSION 


The design analysis of the LP turbine has been carried 
out using the finite elements. A flat triangular shell ele- 
ment and bar elements have been chosen in discretizing 
the casing. The computer code developed uses frontal 
solution procedure since this is very advantageous for 
closed type of structures. The analysis has been done for 
loads which result during lifting of the casing. The com- 
puted results have been compared with experimental 
results and were found to be in good agreement for 
normal lifting. Using this computer code, one particular 
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arrangement of the stiffeners was checked for stresses 
which would result during improper lifting of this struc- 
ture. The stresses induced are within the safe allowable 
limit indicating that the stiffener arrangements are 
adequate. The finite element modeling help evaluate the 
stresses induced in the turbine casing for validating the 
design. Also, the computer code developed can be used 
for the design analysis of box type power plant equip- 
ment such as condensers, large D.C. and A.C. machine 
frames, etc. 
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